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We provide a prescription to construct a microscopic model for heavy lanthanide systems such as Yb 
and Tm compounds by exploiting a j-j coupling scheme. Here we consider a situation with a large spin- 
orbit coupling, in which j=5/2 sextet is fully occupied, while j=l/2 octet is partially occupied, where j 
denotes total angular momentum. We evaluate crystalline electric field potentials and Coulomb interactions 
among the states of the j=7/2 octet to construct a local Hamiltonian in the j-j coupling scheme. Then, it 
is found that the local /-electron states composed of the j=7/2 octet agree quite well with those of seven / 
orbitals even for a realistic value of the spin-orbit coupling. As an example of the application of the present 
model, we discuss low-temperature multipole states of Yb- and Tm-based filled skutterudites by analyzing 
multipole susceptibility of the Anderson model in the j-j coupling scheme with the use of a numerical 
renormalization group technique. ^From the comparison with the numerical results of the seven-orbital 
Anderson model, it is concluded that the multipole state is also well reproduced by the j-j coupling 
model, even when we include the hybridization between conduction and / electrons for the realistic value 
of the spin-orbit coupling. Finally, we briefly discuss future applications of the present prescription for 
theoretical research on heavy lanthanide compounds. 
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1. Introduction 

It has been widely recognized that emergence of heavy 
electron state is understood from the competition in /- 
electron duality nature, 1 ' i.e., itinerancy due to Kondo ef- 
fect 2 ^' vs. localization due to Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction. 5 ~ 7) There appears a quantum 
phase transition at zero temperature in the /-electron state be- 
tween itinerant and localized regions. In such a competing re- 
gion, we have frequently observed unconventional supercon- 
ductivity and non-Fermi liquid behavior due to the effect of 
quantum critical fluctuations. 8 ~ 14) The emphasis on quantum 
critical nature has been summarized in the famous Doniach's 
phase diagram, 15 ' which has been a guiding principle for a 
long time to understand several kinds of anomalous electronic 
properties of heavy fermion materials. 

The concept of quantum critical point seems to be univer- 
sal, since it holds for p-, d-, and /-electron systems. It is an 
important issue to accumulate the experimental facts which 
can be explained by the universal concept of quantum criti- 
cality, although it is sometimes difficult to control experimen- 
tally quantum criticality in actual materials. However, if we 
ignore individual characters of electrons such as orbital de- 
gree of freedom, in general, difference between d and / elec- 
trons cannot be understood only from the control of the ratio 
of Coulomb interaction and electron bandwidth. It seems to be 
also important to emphasize individuality of electron in mate- 
rials, in particular, when we attempt to synthesize new exotic 
and functional materials. 

In this context, we are interested in Yb and Tm systems, 
which have attracted renewed attention due to difference in 
quantum critical nature between Ce and Pr compounds. For 
trivalent ions, one and two / electrons are included in Ce 3+ 
and Pr 3+ , respectively, while one and two / holes exist on 
Yb 3+ and Tm 3+ , respectively. If we use the electron-hole 
symmetry, we expect similar electronic properties between 



/-electron and hole systems. Such a discussion may be also 
found in Pr and Tm compounds. 

Thus far, unconventional superconductivity has been found 
in Ce-based materials since the pioneering discovery of super- 
conductivity in CeCu2Si2- I6) Recently, relatively high super- 
conducting temperature over 2 K has been also observed in 
a material group of CeTIns (T=Ir, Rh, and Co). I7 ~ 19) Then, 
if we simply believe the electron-hole picture, at the first 
glance, it seems to be easy to find superconductivity in Yb 
system. However, in spite of much effort to seek for super- 
conductivity in Yb-based heavy-fermion materials, it has been 
very difficult to synthesize superconducting Yb compounds. 
Recently, superconductivity has been observed in /3-YbAlB4 
with a superconducting temperature T c =80 mK. M) It has 
been claimed that this material exists just on the quantum crit- 
ical point at ambient pressure. In TmsRhgSnig, superconduc- 
tivity has been also found with T C =2.2K. 23) Peculiar reentrant 
properties have been considered to be related to the coexis- 
tence of magnetism and superconductivity. 

If we emphasize similarity between Ce and Yb compounds, 
we prefer to exploit the concept of quantum criticality on the 
basis of the electron-hole symmetry. In this case, different 
quantum critical nature will be explained by the difference in 
the local interactions and the hybridization between conduc- 
tion and localized / electrons. However, one may have a sim- 
ple question whether quantum criticality of one /-hole system 
is really the same as that of one /-electron system. In order 
to clarify /-electron state in Yb compounds, we can choose 
an alternative way to focus on the difference in relevant /- 
electron orbital. This point will be also related to similarity 
and difference between Pr- and Tm-based compounds. 

For the research along such a direction, it is necessary to 
define the /-electron state by a conventional way to include 
many-body effects. Here we recall a couple of schemes for 
the description of local /"-electron configuration, where n 
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denotes the number of / electrons on a localized rare earth or 
actinide ion. One is an LS coupling scheme, in which we con- 
struct the spin S and angular momentum L by following the 
Hund's rules as S=Y^ =i s« and £=X^=i wnere s i an d 
are spin and angular momenta for i-th / electron, respectively. 
As is well known, the Hund's rules are based on the Pauli prin- 
ciple and Coulomb interactions among / electrons. After the 
formation of S and L, the effect of spin-orbit interaction is 
included in the form of S,L-S, where £ is the spin-orbit cou- 
pling in the LS coupling scheme. We note that £>0 for n<7, 
while £<0 for n>7. Note also that a good quantum number 
to label such a state is the total angular momentum J, which 
is given by J=L+S. Then, the ground state is characterized 
by J=\L-S\ for n<7, while J=L+S for n>7. 

As is understood from the above discussion, the LS cou- 
pling scheme is quite useful for the case in which the Hund's 
rule coupling is much larger than the spin-orbit interaction, 
since S and L are formed by the Hund's rule coupling prior 
to the inclusion of the spin-orbit interaction. This assumption 
is considered to be valid for insulating compounds with lo- 
calized / electrons. However, when the spin-orbit interaction 
is not small compared with the Hund's rule coupling, for in- 
stance in actinide compounds, the above assumption is not 
always satisfied. In addition, if the / electrons begin to be 
itinerant due to hybridization with the conduction electrons, 
the effect of Coulomb interactions would thereby be effec- 
tively reduced. In rough estimation, the effective size of the 
Coulomb interaction may be as large as the bandwidth of / 
electrons, leading to a violation of the assumption required 
for the LS coupling scheme. 

For /-electron systems in which the spin-orbit interaction 
becomes larger than the effective Coulomb interactions, we 
prefer to exploit a j-j coupling scheme. 24 ' 25) Here we em- 
phasize that the j-j coupling scheme is convenient for the in- 
clusion of many-body effects by using the standard quantum- 
field theoretical techniques, since individual /-electron states 
are clearly defined, as we explain below. First, we include the 
spin-orbit coupling so as to define the state labelled by the 
total angular momentum j i for the i-th / electron, given by 
ji=Si+li. For / orbitals with 1=3, we obtain an octet with 
j=7/2 and a sextet with j=5/2, which are well separated by 
the spin-orbit interaction. Note that the level for the octet is 
higher than that of the sextet. Then, we consider the effect of 
Coulomb interactions to accommodate n electrons among the 
sextet or octet, leading to the ground state in the j-j coupling 
scheme. For the models of Ce and Pr compounds, the sextet 
should be used for the construction of the effective model. 24 ' 
On the other hand, for Yb and Tm materials, since the sextet is 
fully occupied, we consider the octet to construct the model. 

In this paper, we develop a prescription to construct a mi- 
croscopic effective model for heavy lanthanide systems such 
as Yb and Tm compounds on the basis of the j-j coupling 
scheme. Then, it is shown that the local /-electron state in the 
j-j coupling scheme agrees quite well with that of the origi- 
nal model including seven / orbitals, even for a realistic value 
of the spin-orbit coupling. Next we consider the impurity An- 
derson models to discuss low-temperature multipole states of 
/ electrons. Here we pick up Yb- and Tm-based filled skut- 
terudites as typical examples. The models are analyzed with 
the use of a numerical renormalization group technique. It 
is found that the multipole state is well reproduced by the 



j-j coupling scheme in comparison with that of the seven- 
orbital Anderson model, even when we include hybridization 
between conduction and localized / electrons for the realistic 
value of the spin-orbit coupling. 

The organization of this paper is as follows. In Sec. 2, 
we discuss the local Hamiltonian in the j-j coupling scheme 
in comparison with the results of the original seven-orbital 
model. In Sec. 3, we set the impurity Anderson models for 
filled skutterudites. Then, we show the numerical results for 
multipole susceptibility to discuss the validity of the j-j cou- 
pling model. In Sec. 4, we provide a few comments on future 
issues on the present prescription and summarize this paper. 
Throughout this paper, we use such units as fc B =/j=l. 

2. Local Hamiltonian in a j-j Coupling Scheme 

2.1 Original Seven-Orbital Model 

In general, the local /-electron Hamiltonian is given by 

H\oc = H so + HqEF + H int . (1) 

The first term denotes the spin-orbit coupling, given by 

H S o A ^ ' Cm.cr-.m' ,a' fmcrfm'cr' i (2) 
m.(7,m' .a' 

where A is the spin-orbit interaction, f ma is the annihilation 
operator of / electron, a=+l (—1) for up (down) spin, m is 
the z-component of angular momentum £=3, and the matrix 
element ( is given by 

Cm.a;m.a — 77l(j/2, 

( m +a,-o; m ,a = y/ l{l + 1) - m{m + <j)/2, 

and zero for other cases. 

The second term denotes crystalline electric field (CEF) po- 
tential, given by 

-HCEF = ^ B m>m > fmvfm'a, (4) 

where B m ^ m i is determined from the CEF table for 
J=£=3. 26 ' 27 1 Note that electrostatic CEF potentials do not act 
on /-electron spin. Since we will consider later the multipole 
state of heavy lanthanide filled skutterudites, here we show 
Bm,m' of the cubic system with Ti, symmetry. 28) The results 
are given by 

£?3, 3 = £_3,_3 = 180B 4 + 180S 6 °, 

B 2 ,2 = £-2.-2 = -420B2 - 10805°, 

B hl = S_i _j = 605^ + 2700Bg, 

B 0fl = 360B2 - 360073g, 

B 3! _ 1 =B_ 3tl = 12VT5(Bi + 5Bl), 

B 2 - 2 = 6051 - 360B|, K ' 

B 3 ,i = £-3,-i = 24V15BI, 

B 2 .o = £-2,0 = -48V30BI, 

B lt -! = 360B|, 

B 3 !-3 - 360BI 

Note the relation of B m ^ m i=B m ^ m . We also note the rela- 
tions of B\=hB%, B\=-2\B%, and B%=-B%. Following the 
traditional notation, 29 ' we define 

Bl = Wx/F(4), 

B$ = W(l-\x\)/F(6), (6) 
Bl = Wy/F\Q), 

where x, y, and the sign of W specify the CEF scheme 
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for Th point group, 28 ' while the absolute value of W deter- 
mines the energy scale of the CEF potential. Concerning non- 
dimensional parameters, 5(4) and 5(6), we choose F(4)=15, 
F(6)=180, and 5*(6)=24 for J=3. 

Finally, 5; nt denotes Coulomb interaction term, given by 

Hi n t= ^ ' ^ ' -fmim2 ! m 3 m4/m|j/m, ff '/m3g'/m4g; (7) 

where the Coulomb integral l m , m!;m3m4 is expressed by 

6 



1 mim2 ; m3m4 



k=0 



F k c k {m 1 , m 4 )c fe (TO 2 , m 3 ). 



(8) 



Here 5 fc is the Slater-Condon parameter 30 ' 31 ' and Cfc is the 
Gaunt coefficient 32,33 ' which is tabulated in the standard text- 
books of quantum mechanics. 34 ' Note that the sum is limited 
by the Wigner-Eckart theorem to k=0, 2, 4, and 6. 

2. 2 Effective Hamiltonian for j= 7/2 Octet 

In order to obtain the model in the j-j coupling scheme, 
we transform the /-electron basis between (to, a) and (j, ^) 
representations, connected by Clebsch-Gordan coefficients, 
where j is the total angular momentum and /i is the z- 
component of j. When we define as the annihilation op- 
erator for / electron labelled by j and [i, the transformation is 
given by 

m,<7 

where the Clebsch-Gordan coefficient Cj.^ [m . a is give by 



C, 



5/2,fi;fi— o j 2,a 



-oV (7/2 - a^)/7, 



C 7 /2, W([1 - CT /2,a = ^(7/2 + <7//)/7, 



(10) 



and other components are zero. 

After the transformation, the spin-orbit coupling term is di- 
agonalized as 



(11) 



sider only the j=7/2 octet. 

The local model in the j-j coupling scheme is given by 



H 



loc 



He 



EF 



Hi, 



(14) 



where Hqef is the CEF potential in the j=7/2 octet, given as 
HcEF = Y,B ll>v flf v . (15) 

Here is the creation operator of / electron in the ^i-state 
and n indicates the z-component of the total angular momen- 
tum which specifies the state in the j=7/2 octet. The CEF po- 
tential in the j=7/2 octet is given by 27 ' 



7/2,7/2 
5/2,5/2 
3/2,3/2 



#1/2,1/2 " 
-87/2,-1/2 



5_ 
B_ 
B_ 
B 



7/2,-7/2 
5/2,-5/2 
3/2,-3/2 
-1/2,-1/2 



= 5_ 



7/2,1/2 



-65/2,-3/2 = -B-5/2,3/2 = 60^3(5, 



42052 + 12605g, 
-7805° - 63005g, 
-18052 + 113405°, 
54052 - 63005§, 

f 155 4 ), 
75|), 



12V35(5 4 



(16) 



5. 



7/2,3/2 



5_ 



7/2,-3/2 



120V215|, 



#5/2,1/2 — £-5/2,-1/2 — -504\/55g, 



5. 



3/2,-1/2 



5_ 



3/2,1/2 



168V155 2 , 



-87/2,-5/2 = £-7/2,5/2 = 360v / 7"-E?6> 

where we note again the relations of 5 4 =552, 5 4 =— 215g, 
and Bq=—Bq. The CEF parameters for j=7/2 are related to 
those for J=£=3 as 



£ 4 ° = (/W&)#4° = 35 4 °/7, 
£ 6 ° = (77/2/73)5 6 - 5 6 °/7, 
£1 = (77/2/73)5 6 2 = 5 2 /7, 



(17) 



where /3j and 77 are fourth- and sixth-order Stevens fac- 
tors, respectively. 26 ' We note that /3 3 =2/495, /5 7/2 =2/1155, 
73 =-4/3861, and 77/2 =-4/27027. 

The second term in eq. (14) indicates the Coulomb interac- 
tions in the j=7/2 octet, which is given by 



Hn 



E 

1-1,1;. fi f 



1 



(18) 



with A 5 /2=— 2A and A 7 /2=(3/2)A. The CEF and Coulomb in- 
teraction terms are, respectively, given by 



(12) 



Jl^l,j2M2 



and 



Jl~j4 Ml~M4 



where 5 and / are the CEF potential and Coulomb interac- 
tion, respectively, in the basis of j and fi. 

In the present paper, we consider the model for heavy lan- 
thanide systems with n>7. In the limit of large A for the j-j 
coupling scheme, j=5/2 sextet is fully occupied, while j=7/2 
octet is partially occupied. Thus, here we simply discard all 
the j=5/2 states and keep only the j=7/2 octet. Namely, we 
accommodate n— 6 electrons in the j=l/2 octet. Note that in 
this approximation, the spin-orbit coupling is given by the ef- 
fect of potential energy which does not depend on the orbitals. 
Since such an energy can be included in the chemical poten- 
tial shift, we do not consider explicitly 5 so in the following. 
Hereafter, we suppress the subscription j in since we con- 



where / is the matrix element for Coulomb interactions 
among 7=7/2 states. In order to classify the Coulomb inter- 
actions in the j=7/2 octet, we consider the situation where 
we accommodate two electrons in the octet. Note that the al- 
lowed values for total angular momentum J are 0, 2, 4, and 
6 due to the Pauli principle. Thus, the Coulomb interaction 
term should be written in a 28x28 matrix form. Note that 
"28" is the sum of the basis numbers for singlet ( J=0), quin- 
tet (J=2), nonet (J=4), and tridectet (J=6). As is easily un- 
derstood, this 28 x28 matrix can be decomposed into a block- 
diagonalized form labelled by J z , including one 4x4 matrix 
for J z =0, four 3x3 matrices for J z =±l and ±2, four 2x2 
matrices for J z =±3 and ±4, and four lxl for J z =±5 and 
±6. We skip the details of tedious calculations for the evalu- 
ation of matrix elements and show only the results by using 
the parameters Ek (fc=0,l,2,3), 35 ' 36 ' which are related to the 
Slater-Condon parameters F k as 37 ' 



E a = F° 25-^2 5_^. 

1? _ 40. F 2 5 l 7 _8_ P4 2 ! 1 200 p6 
1 567 231 T H583 

5o = -^-F 2 - -2-5 4 

1617 5929 1 > 

E 3 = J^F 2 + J^*F A 



125 



J. i:>n:s 



(19) 



100 



43659 



17787 



1656369 



F b . 
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For J z =6 and 5, we obtain 

4/2,5/2;5/2,7/2 = S - 154S 3 , 

and 

4/2,3/2;3/2,7/2 = E — 154S 3 , 

respectively. For J Z =A and 3, we obtain 



'J/2,l/2;l/2,7/2 
4/2,3/2;3/2,5/2 
4/2,l/2;3/2,5/2 



So - 35S 2 /2 - 119S 3 /2, 
So - 75£ 2 /2 + 97S 3 /2, 
Vl05(5S 2 - 27S 3 )/2, 



and 



'7/2,-l/2;-l/2,7/2 



£ - 35S 2 + 35S 3 



4/2,l/2;l/2,5/2 — S - 20S 2 - 46S 3 , 
4/2,-l/2;l/2,5/2 = \/7(10S 2 - 54S 3 ). 



For J z =2 and 1, we obtain 



'7/2,-3/2;-3/2,7/2 



'5/2,-l/2;-l/2,5/2 



So 

So 



21S 2 
35S 2 



98S 3 , 
46S 3 , 



4/2,1/2:1/2.3/2 — S + 30S 2 + 8OS3, 

= -\ZT05(2S 2 +9S 3 ), 

9x/35(S 2 - S 3 ), 
-V3(25S 2 +63S 3 ), 



4/2,-3/2;-l/2,5/2 
4/2,-3/2;l/2,3/2 = 
4/2,-l/2;l/2,3/2 " 



and 



'_7/2,-5/2;-5/2,7/2 
4/2,-3/2;-3/2,5/2 



3/2,-1/2; 
4/2,-5/2; 



-1/2,3/2 
-3/2,5/2 



S + 49S 2 /2 + 259S 3 /2, 

S + 27S 2 + 17S 3 , 

S - 15S 2 /2 - 29S 3 /2, 

-V2T(13S 2 + 9S 3 ), 



4/2,-5/2;-l/2,3/2 — 9V105(S 2 - S 3 )/2, 
4/2,-3/2;-l/2,3/2 = ~3\/5(S 2 +21S 3 ), 

Finally, for J z =0, we obtain 



4/2,- 


-7/2 


-7/2,7/2 


= S + Si + 49S 2 


+ 105S 3 , 


4/2,- 


-3/2 


-3/2,3/2 


= So H- Si — 29S2 


+ 51S3, 


4/2,- 


-5/2 


-5/2,5/2 


= So + Si + 9S2 — 


75S 3 , 


4/2,- 


-1/2 


-1/2,1/2 


= S + Si + 15S 2 


+ 51S3, 


4/2,- 


-7/2 


-5/2,5/2 


= -Si - 49S 2 /2 + 


49S 3 /2, 


4/2,- 


-7/2 


-3/2,3/2 


= Si - 21S 2 - 56S 3 , 


4/2,- 


-7/2 


-1/2,1/2 


= -Si + 105S 2 /2 - 


f 49S 3 /2, 


4/2,- 


-5/2 


-3/2,3/2 


= -Si + 21S2/2 - 


I3IS3/2, 


4/2,- 


-5/2 


-1/2,1/2 


= Si + 15S 2 /2 - 92S 3 , 


4/2,- 


-3/2 


-1/2,1/2 


= -Si - 45S 2 /2 - 


131S3/2, 


Note here the following relations: 





*/^,^;^',/y 



and 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



(28) 



By using these two relations and eqs. (20)-(26), we can obtain 
all the Coulomb matrix elements. 

2. 3 CEF Energy Levels 

Let us now consider the situation in which two electrons 
are accommodated in the j=H2 octet. This situation indi- 
cates the case with 8 electrons in / orbitals, corresponding 
to / 8 configuration of Tb 3+ ion. When we diagonalize the 
28 x 28 matrix for Coulomb interaction terms, we can easily 
obtain the eigen energies as S — 154S 3 for the J=6 tridectet, 



S - 55S 2 + 143S 3 for the J=4 nonet, S + 99S 2 + 143S 3 
for the J=2 quintet, and S + 4Si for the J=0 singlet. These 
values are exactly the same as those obtained in the nuclear 
shell theory in the j-j coupling scheme. 38 ' For typical values 
of Slater-Condon parameters, we find that the ground state is 
specified by J=6 in the j-j coupling scheme. For Tb 3+ ion, in 
the LS coupling scheme, we obtain the ground-state level as 
7 F with 5=3 and L=3 from the Hund's rules. On further in- 
clusion of the spin-orbit interaction, the ground state becomes 
characterized by J=6, expressed as 7 Fq in the traditional no- 
tation. Note that we are considering a two-electron problem. 
Thus, when we correctly include the effects of Coulomb inter- 
actions, the same quantum number as that in the LS coupling 
scheme is obtained in the j-j coupling scheme for the ground- 
state multiplet. 

In order to discuss the CEF energy levels, it is necessary to 
determine the values of local interactions. Among them, con- 
cerning the Slater-Condon parameters, we set S°=10 eV by 
hand. The magnitude of S° is related to the absolute value 
of the ground state energy. It can be evaluated by the first- 
principles calculation, but it is out of the scope of the present 
paper. Other Slater-Condon parameters are determined so as 
to reproduce excitation spectra of Pr 3+ ion. 39 ' 40) After the 
fitting, we obtain S 2 =8.75 eV, S 4 =6.60 eV, and S 6 =4.44 
e y.4i,42) As long as we ignore the difference in lanthanide 
ions, e.g., the size of ion radius, we use these values for all 
lanthanide ions. On the other hand, as for the spin-orbit cou- 
pling A, we use the value which has been determined experi- 
mentally for each lanthanide ion. 43 ' 

Concerning CEF parameters, it is necessary to specify the 
actual material, since they depend on the crystal structure and 
the kinds of ligand ions. Here we consider the case of filled 
skutterudite materials, since we will analyze the multipole 
state of filled skutterudites later in this paper. The CEF pa- 
rameters are different from material to material even if we fix 
the material group, but the typical values are W=— 0.4 meV, 
y=0.3, and .t=0.3, which are determined so as to reproduce 
quasi-quartet CEF scheme of PrOs 4 Sbi 2 . 44 "* 6) Note that for 
Pr atom, we use A=0.095 eV from the experimental value. 

In Fig. 1(a), we show the results of CEF energy levels vs. 
x for the case of n=2, corresponding to Pr 3+ , with the use of 
the above parameters. The vertical dash line denotes the po- 
sition of £=0.3 and we can understand that the ground state 
is singlet and the first excited state is T^ 2 ^ triplet with 
the small excitation energy. This is considered to be a typi- 
cal situation of PrOs4Sbi 2 . When we change the values of x 
and/or y, we can obtain another situation for different filled 
skutterudite material. 

In Fig. 1(b), we show the CEF energy levels vs. x for the 
case of n=13 with W=— 0.4 meV and y=0.3. The spin-orbit 
coupling A is set as 0.356 eV for Yb atom. 43) In the O h point 
group, the octet is known to split into two doublets and one 
quartet. 29 ' In the Xh point group, on the other hand, two dou- 
blets are mixed. 28 ' In the present parameters, we always find 
the doublet ground state, irrespective of the values of x. 

In Figs. 2, we show the results for the case of n=12 cor- 
responding to Tm 3+ ion, in order to see the validity of the 
j-j coupling scheme. In Fig. 2(a), the result for the original 
seven-orbital model in eq. (1) is shown. In the Oh group, the 
tridectet of J=6 is split into two singlets, one doublet, and 




Fig. 1. CEF energy levels for (a) n=2 and (b) n=13. The vertical line de- Fig. 2. CEF energy levels vs. x for n=12 (a) in the original seven-orbital 
note the position of x=03. Concerning other parameters, see the main text. model eq. (1) and (b) in the j-j coupling model eq. (14). 



three triplets. ' In the Xh group, on the other hand, two sin- 
glets are mixed. Three triplets are also mixed due to the y- 
term. In the present parameters, the ground state around at 
£=0.3 is characterized by singlet. However, when we in- 
crease the value of x, we find the change of the ground state 
from singlet to rjj non-Kramers doublets. Such a change 
is found to occur around at x ~ 0.9. 

In Fig. 2(b), we show the CEF energies vs. x of the j-j 
coupling model eq. (14) with the use of the same parame- 
ters, except for the value of the spin-orbit coupling. In the j-j 
coupling scheme, A is set as infinity. Nevertheless, in the first 
impression, even if the actual value of A is finite, the results of 
the j-j coupling model agree quite well with those of the orig- 
inal seven-orbital model in Fig. 2(a). If we effectively change 
the absolute value of W as a fitting parameter, Fig. 2(a) can 
be reproduced by the j-j coupling model quantitatively. 

The reason why the CEF energy levels are well reproduced 
by the j-j coupling scheme even for the finite value of A 
is as follows. Since the maximum value of the difference in 
the z-component of total angular momentum is seven among 
j=7/2 states, the sixth-order CEF potential can be included 
in the j-j coupling scheme, in sharp contrast to the case of 
j=5/2. 24) Thus, the difference in the value of A does not pro- 
vide serious effect on the CEF ground state, as long as we 
consider A^>| W\. In fact, after lengthy algebraic calculations, 
we obtain the overlap integral of the ground states for n=12 
as ($|$ J _ J )= v /677=0.926, 35 - 36) where |$) and de- 



note the CEF ground state of the seven-orbital model and the 
j-j coupling one, respectively. 

We do not show further the results for the cases of n<12, 
but the CEF states of the j-j coupling model can reproduce 
well those of the original seven-orbital model. Thus, we con- 
clude that the local /-electron state of heavy lanthanide sys- 
tems is well approximated by the j-j coupling scheme. This 
is one of important messages of the present paper. 

3. Numerical Analysis of Impurity Anderson Models 

We have explained the prescription to obtain the local 
Hamiltonian in the j-j coupling scheme. Even if we use only 
the j=ll2 octet, it is possible to reproduce local multi-/- 
electron state, which agrees well with those obtained in the 
original seven-orbital model. By including further the itiner- 
ancy of / electrons, we can discuss magnetism and supercon- 
ductivity of heavy lanthanide compounds from a microscopic 
viewpoint. In this section, as an example, we discuss the mul- 
tipole state of Yb- and Tm-based filled skutterudites by using 
the Anderson models. Then, we show the effectiveness of the 
j-j coupling model in the microscopic level. 

3.1 Anderson Models 

The seven-orbital Anderson Hamiltonian is given by 

H=J2 £ k c L c ka+ (^™ c L/«"x + h - c -)+#ioc, (29) 
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where denotes conduction electron dispersion, Cfc CT indi- 
cates the annihilation operator for conduction electron with 
momentum k and spin a, V m is the hybridization between 
conduction and / electrons, and the local /-electron term 
H\ oc is already given in eq. (1). For filled skutterudites, 
the main conduction band is given by a u , constructed from 
p-orbitals of pnictogen. 47 ' Note that the hybridization oc- 
curs between the states with the same symmetry. Since the 
a u conduction band has xyz symmetry, we set V 2 =V/\/2, 
V_ 2 =-V/V2, and zeros for other values of m. The hybridiza- 
tion is fixed as V=0.05 eV and a half of the bandwidth of a u 
conduction band is set as 1 eV. 

The j-j coupling Anderson model is given by 

H = Y, £ kc{ a c ka + (V^AaU + h - c -) + Hioc, (30) 



where the local /-electron term in the j-j coupling scheme 
H\ oc is given by eq. (14). Since the a u conduction band has 
xyz symmetry, which is described by in the T h group, i.e., 
Tj in the O h group, we set V^^/ 2 =V\/3/2, V^._ 3 / 2 = — 
V\.,-5/2=— Vy/3/2, Vj.. 3/2=^/2, and zeros for other cases. 
Note that the connectivity between /i and a is determined by 
the definition of pseudo-spin of /-electron state on the basis 
of the time reversal symmetry. For the j-j coupling model, 
we also set F=0.05 eV for simplicity, since this value is not 
so important at low enough temperatures in the following dis- 
cussion. Note that for the Th group for filled skutterudites, 
and T-T doublets in the Oh group are mixed and they are 
expressed as T 5 ^ and T 5 t " 2 \ 

It should be noted that we do not show explicitly the chem- 
ical potential terms both in the models, but in actual calcula- 
tions, we set the value of the chemical potential so as to fix 
the local /-electron number as n=13 or 12. 

3.2 Multipole Operator 

In order to discuss the multipole state, it is necessary to 
define the multipole operator. The details can be found in 
Refs. 41 and 42, but here we briefly explain the method to de- 
fine the multipole operator to make this paper self-contained 
with some additional comments. 

When we consider multipole operator for / electrons, it 
should be defined in the one-body form as an extension of 
charge and total angular momentum operators on the basis of 
a belief that the multipole denotes the combined degree of 
freedom of spin and orbital 42) 

In the multipole expansion of potential in electromag- 
netism, higher electric and magnetic multipole moments ap- 
pear in the coefficients of the expansion by the spherical har- 
monics Ylm with larger angular momentum. In group theory, 
Ylm is defined by the basis of irreducible representation 
of the rotation group R, expressed as 



RY LM =J2 Ylm ' D 



(L) 
MM' 



(31) 



M' 



In order to define /-electron multipole operator, on the anal- 
ogy of the multipole expansion, we exploit a concept of spher- 
ical tensor operator in the quantum mechanics of angular mo- 
mentum. 48 ' When we consider the rotation of operator T, 
we obtain a set of operators f^={f^ k) } with (2k + 1)- 



components (q = —k, — k + 1, • • • , k — 1, k), given by 

RT^R- 1 =Y ?P D ( q k J . (32) 



~ (k) 

Namely, Tq is transformed like a basis of irreducible repre- 
sentation for the rotation. Such Tq^ is called spherical 
tensor operator of rank k. 

Thus far, we have implicitly assumed /-electron density 
in an isolated ion, but in actuality, rare-earth ions are put in 
the crystal structure. Then, it is convenient to change from 
spherical to cubic tensor operators, given by 



(33) 



where k is a rank of multipole, an integer q runs between — k 
and fc, 7 is a label to express Oh irreducible representation, 
and G^} q is the transformation matrix between spherical and 
cubic harmonics. Then, the cubic tensor operator for / elec- 
tron is expressed in the second-quantized form as 



f (fc) = \p T (fc, 7 ) f t f 
7 / j mcr,m'u' J maJn 



(34) 



raa.nv a' 



Throughout this paper, we use the cubic tensor operator as 
multipole. 

As for the classification of multipole, we use the notations 
in the group theory. We express the irreducible representation 
of the CEF state by Bethe notation in this paper, but for multi- 
poles, we use short-hand notations by the combination of the 
number of irreducible representation and the parity of time re- 
versal symmetry, g for gerade and u for ungerade. Note also 
that for the Th group, Y\ and T 2 of Oh are mixed. We remark 
that T4 and T 5 of Oh are also mixed in Th. Thus, we obtain 
six independent multipole components as lg+2g, 2u, 3g, 3u, 
4g+5g, and 4u+5u for filled skutterudites. Note that lu does 
not appear within rank 7. 

The coefficient T,i fc ' 7 2/ / is calculated from the spherical 
tensor operator as follows. First we change the /-electron ba- 
sis from (m, a) to (j, p). Note that j takes 7/2 and 5/2 for / 
electrons. For a certain value of angular momentum j and its 
z-component /j,, the matrix element of spherical tensor opera- 
tor is easily calculated by the Wigner-Eckart theorem as 

(jp\T^\jp') = OI )f fc _! ll 1 j) 0>b> / fcg), (35) 
V2j + 1 

where (J M\J' M' J" M") denotes the Clebsch-Gordan co- 
efficient and (j||T( fc )||j) is the reduced matrix element for 
spherical tensor operator, given by 



(j\\T^\\j) = 



_1_ / (2j + fc + 1)! 
2 fc \/ (2j-k)\ ' 



(36) 



Note that k < 2j and the highest rank is 2j. The coefficient 
Tmfm'a' ' s obtained by returning to the basis of (m, a) from 
(j, p). The final result is given by 



r (fe,7) 

-ma,m'cr' / t ^7i<? ^/2j + 1 



(37) 



x (jp\tms'^)(jp.'\tm's C ^), 
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where 1=3, s=l/2, j=£±s, and fj, runs between —j and j. 

For the j-j coupling scheme in the j=7/2 octet, we should 
discard the contribution from j=5/2 sextet. Then, the multi- 
pole operator in the j-j coupling scheme is expressed in the 
second-quantized form as 



5.0 



f(k) _ f,(fe,7) ft f , 



(38) 



A*, A*' 



where // denotes the z-component of j=7/2. The coefficient 



T. 



01|r (fc) lb 
V2JTT 



'-{jlAw'H, (39) 



where j is fixed as j=7/2 in this equation. We use this defini- 
tion for the calculation of the multipole susceptibility in the 
j-j coupling scheme. 

It should be noted here that multipoles belonging to the 
same symmetry are mixed in general, even if the rank is dif- 
ferent. In addition, multipoles are also mixed due to the effect 
of CEF potentials of the group. Namely, the /-electron 
spin-charge density should be given by the appropriate super- 
position of multipoles, expressed as 



(40) 



„(*) 



In order to determine the coefficient p\ ' , it is necessary to 
evaluate the multipole susceptibility in the linear response 
theory. However, multipoles belonging to the same symme- 
try are mixed in general, even if the rank is different. In addi- 
tion, multipoles are also mixed due to the CEF effect. Thus, 

(k) 

it is natural to define p\ ' by the eigenstate of susceptibility 
matrix 

1 p-Et/T _ p-Ej/T 

Xkl ^, =i y - l (ijprw - 



z ^ 

hi 

x <7|[Tf > 



Ei 



id 



(41) 



where Ei is the eigenenergy for the i-th eigenstate \i) of H 
or H, T is a temperature, e~ E * /T (i\f( k) \i)/Z, and 

Z is the partition function given by Z=£\ e- Ei / T . Note that 
the multipole susceptibility is given by the eigenvalue of the 
susceptibility matrix. 

3.3 Method 

In order to evaluate the multipole susceptibility of the im- 
purity Anderson model, here we employ a numerical renor- 
malization group (NRG) method. 49 ' In this technique, we can 
include efficiently the conduction electrons states near the 
Fermi energy by discretizing momentum space logarithmi- 
cally. Note that in actual calculations, it is necessary to in- 
troduce a cut-off A for the logarithmic discretization of the 
conduction band. Due to the limitation of computer resources, 
we keep only M low-energy states. In this paper, we set A=6 
and M=2000. The temperature T is defined as T=A~( Ar - 1 )/ 2 
in the NRG calculation, where N is the number of the renor- 
malization step. With the use of NRG technique, we evaluate 
entropy <Si mp , specific heat C\ mp , and multipole susceptibil- 
ity x- In particular, the optimized multipole state is defined 



1.0 



I 1 1 1 lll| — I I 1 1 1 lll| — I I 1 1 1 ■ 1 1 1 — I I 1 1 1 ■ 1 1 1 — I I 1 1 1 lll| — I I 1 1 1 lll| — I I 1 1 1 ■ 1 1 1 — I I 1 1 1 ■ 1 1 1 — I I 1 1 Mii| — I I I n 

(a) original 

■— ■ entropy 
* specific heat 



4.0 


_ n=13 


CD 


x=0.3 


J 3.0 


" y=0.3 

_ W=-0.4meV 


and 




J°2.0 






a 

E 




00 





0.0 
10 

5.0 




10 io 9 io 8 id 7 io 6 io 5 10 4 16 3 io 2 io 1 10 

T (unit of eV) 

I I iinij — I 1 1 1 1 ni| — I 1 1 1 1 ni| — i i i ■ ■ ■ ■ 1 1 — i i i ■ ■ ■ ■ 1 1 



4.0- 



Q E 3.0 

CO 



2.0 



E 
00 



1.0 



0.0 



l| — I I 1 1 1 lll| — I I 1 1 1 ■ 1 1 1 — I I 1 1 1 ■ 1 1 1 



(b) j-j 

n=13 
x=0.3 
y=0.3 

W=-0.4meV 



entropy 
specific heat 



log 8 



log 2 




i i i i i ' ' 



_-9 .„-S tZH ,„-6 ,„-5 -4 . -3 .„-2 ,„0 

10 10 10 10 10 10 10 10 10 10 10 
T (unit of eV) 

Fig. 3. (Color online) Entropy 5i mp and specific heat Ci m p of (a) the orig- 
inal seven-orbital model eq. (29) and (b) the j-j coupling model eq. (30) 
for n=13 and £=0.3. 



by the eigen state with the maximum eigen value of multipole 
susceptibility matrix eq. (41). 

3.4 Results 

Let us now show our numerical results. First we consider 
the case of n=13. In Figs. 3, we depict the results of entropy 
and specific heat both for the original seven-orbital Anderson 
model eq. (29) and the j-j coupling Anderson model eq. (30). 
Except for the high-temperature region such as 1 eV, we 
do not find significant difference between two panels. Since 
the local /-electron number is fluctuating due to the effect 
of hybridization, the /-electron wave function of the original 
seven-orbital Anderson model is not equal to that of the j-j 
coupling Anderson model. However, as naively expected from 
the similarity in the local /-electron states, the j-j coupling 
model can reproduce well the results of the original model. 

The difference between both models can be observed at 
high temperatures as 1 in the region without enough renor- 
malization steps. Since in the j-j coupling model, we discard 
the j=5/2 states, it is natural that there appears difference from 
the original seven-orbital model at high temperatures. Here 
the energy unit is a half of the conduction bandwidth, which 
is in the order of eV. Since A is in the order of 0. 1 eV, it is 
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Fig. 4. (Color online) Multipole susceptibilities of (a) the original seven- 
orbital model eq. (29) and (b) the j-j coupling model eq. (30) for n=13 
and £=0.3. 



reasonable that the deviation can be found even at T>0.1. 

Around at a temperature of 0.1, we find a plateau in the 
entropy with the value of log 8. It is easily understood that 
this is due to 8-fold degeneracy of j=7/2, since seven elec- 
trons (or one hole) are included in the j=7/2 octet. A part of 
the entropy log 8 is released around at T=0. 01~0. 001 and a 
peak in the specific heat is found at T=0.005. Then, we find 
residual entropy of log 2 in the temperature region less than 
T=0.001. The appearance of the partial entropy release is un- 
derstood due to the CEF energy splitting. Namely, as observed 
in Fig. 1(b), the octet is found to be split in the CEF energy 
scale in the order of 0.01 eV. 

In the present Anderson models for filled skutterudites, we 
consider only the a u conduction band composed of p elec- 
trons in pnictogens. As shown in the explanation of the mod- 
els, this conduction band is hybridized with state in the 
Oh group. In the Th group, Tg and Ty are mixed to form 
two doublets. Namely, the Tj component of the CEF ground 
state for n=13 in the Th group is hybridized with the con- 
duction band, while the component remains even at low 
temperatures. Thus, there occurs residual entropy of log 2 in 
the present model. Of course, when we consider further the 
conduction bands other than a u , the residual entropy should 



rank k 


7 


Seven-orbital model 


j-j coupling model 


1 


4u 


-0.25091 


-0.29008 


3 


4u 


0.18220 


0.18893 


3 


5u 


0.24455 


0.26034 


5 


4u(l) 


0.66679 


0.65432 


5 


4u (2) 


0.11779 


0.11847 


5 


5u 


-0.11770 


-0.11974 


7 


4u(l) 


0.03571 


0.03743 


7 


4u (2) 


-0.57753 


-0.56282 


7 


5u (2) 


-0.19299 


-0.19418 



(k) 

Table I. Coefficients pi, of low-temperature multipole state eq. (40) with 
the largest eigenvalue in the original seven-orbital model eq. (29) and the 
j-j coupling model eq. (30) for n=13 and x=03. 



be finally released at low temperatures. 

Let us move on to the numerical results for multipole sus- 
ceptibility. In Figs. 4, we show the temperature dependence of 
eigenvalues \ of the susceptibility matrix which are classified 
by the symmetry. Note that we plot T\, not x> which is the 
Curie constant for the multipole susceptibility. In this case, 
the difference between both models is not so significant even 
at high temperatures. At low temperatures, the magnitude of 
multipole susceptibility in Fig. 4(a) is slightly different from 
that in Fig. 4(b), but its difference is very small. 

In order to confirm the effectiveness of the j-j coupling 
model, let us turn our attention to the multipole states, not 
the eigenvalues, of the multipole susceptibility. In Table I, we 

(k) 

explicitly list the numbers of the component p\ of the mul- 
tipole state at low enough temperatures for both models. First 
we note that 4u and 5u are mixed due to the effect of the Th 
group and higher-order multipoles are also included with sig- 
nificant weights. In general, there is no explicit relation be- 
tween admixture and rank in the multipole state. It is not sur- 
prising to obtain significant components of higher-order mul- 
tipoles. When we compare the value of each component, of 
course, there exists difference between two models, but we 
can conclude in a satisfactory level that the multipole state 
of the original seven-orbital Anderson model is reproduced 
by the j-j coupling model. Thus, the j-j coupling model is 
useful to analyze the /-electron state with the use of small 
numbers of relevant / orbitals. 

Next we consider the case of n=l2. First we set x as x=Q3 
which is considered to be an appropriate value for filled skut- 
terudites ROs4Sbi2, even if rare-earth atom R is substituted. 
For n=12, as observed in Figs. 2, the CEF ground state is 
singlet. In Figs. 5, we show the numerical results of entropy 
and specific heat both for the original seven-orbital and the j- 
j coupling models. Again we see that both panels agrees well 
with each other, except for difference in the high-temperature 
region T~l. Thus, we reconfirm that the j-j coupling model 
works well also for the case of n=12 with two / holes. 

For a temperature between 0.01 and 0.1, we observe a short 
plateau of log 13 in the entropy, which is considered to be due 
to 13-fold degeneracy of the J=6 state. Then, the entropy of 
log 13 is released to arrive at the singlet ground state. In the 
singlet ground state, we expect no multipole moment. In fact, 
as shown in Figs. 6, multipole susceptibilities vanish at a tem- 
perature at which the specific heat shows a peak due to the re- 
lease of entropy log 13. Note that for the j-j coupling model, 
at high temperatures, we see significant difference in multi- 
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Fig. 5. (Color online) Entropy 5i mp and specific heat Ci mp of (a) the orig- 
inal model seven-orbital eq. (29) and (b) the j-j coupling model eq. (30) 
for n=l2 and £=0.3. 



Fig. 6. (Color online) Multipole susceptibilities of (a) the original seven- 
orbital model eq. (29) and (b) the j-j coupling model eq. (30) for n=12 
and x=0.3. 



pole susceptibilities between Figs. 6(a) and 6(b). Probably it 
is due to the difference in high-energy local states between 
the original seven-orbital and the j-j coupling models in the 
combination with the lack of the renormalization steps. In any 
case, when we further make the renormalization process, we 
finally obtain the same behavior in multipole susceptibility. 

Even at low enough temperatures, no multipole susceptibil- 
ity is observed for n=12 and x=Q3, but it is interesting to con- 
sider a possibility of heavy-electron state in Tm-based filled 
skutterudites. Namely, on the basis of the present numerical 
calculations, we observe large entropy release such as log 13 
at relatively high temperature. It may be risky to conclude 
the heavy-electron state only from the present results, but it 
seems to be interesting to perform the measurements of basic 
bulk properties of Tm-based filled skutterudites, although it 
may be difficult to synthesize actually Tm-based filled skut- 
terudite compounds. 

Let us again turn our attention to multipole state for n=l2. 
Here we increase the value of x by assuming that the value 
of x is controlled experimentally due to the substitution of 
transition metal atoms and/or pnictogens. In Figs. 7 and 8, we 
show the numerical results for n=12 and £=1.0. Note that at 
x=\ .0, the local CEF ground state is r^ 3 non-Kramers doublet 



with the first excited state of triplet, as observed in Figs. 2. 

In Fig. 7(a), we show the results of entropy and specific 
heat. Except for the high-temperature region, we again ob- 
serve that both panels are similar to each other. In this case, 
after a short plateau of log 13 around at T~0.1, we observe 
the remnant of plateau of log 5 due to 5-fold degeneracy of 
quasi-quintet composed of r^ 3 doublet and T± triplet. Then, 
we arrive at the residual doublet state composed of a couple of 
electrons in Tg 7 quartet state (T^ in the G\ group). Since in 
the present model, we consider only the single a u conduction 
band, which is hybridized with iy component of doublet 
state, / electrons in T^ 7 states are considered to be localized. 
As mentioned above, in actual materials, the residual log 2 
entropy should be finally released, since there exist other con- 
duction bands such as e u which hybridize with Tq 7 states. 

Let us explain the results for multipole susceptibility. It is 
observed that except for the high-temperature region larger 
than T=0.1, Figs. 8(a) and 8(b) agree well with each other. At 
high temperatures, significant difference can be found in mul- 
tipole susceptibilities between Figs. 8(a) and 8(b), but it is due 
to the same reasons as those in Figs. 6. At low enough temper- 
atures, we find two kinds of residual multipole states which 
are expected to be dominant in actual materials, although or- 
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Fig. 7. (Color online) Entropy 5i mp and specific heat Ci mp of (a) the orig- 
inal model seven-orbital eq. (29) and (b) the j-j coupling model eq. (30) 
for n=l2 and £=1.0. 



Fig. 8. (Color online) Multipole susceptibilities of (a) the original seven- 
orbital model eq. (29) and (b) the j-j coupling model eq. (30) for n=12 
and x=1.0. 
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Table II. Coefficients pi,, of low-temperature multipole state eq. (40) with 
the largest eigenvalue in the original seven-orbital model eq. (29) and the 
j-j coupling model eq. (30) for n=12 and x=l.0. 



(k) 

Table III. Coefficients p~ f of low-temperature multipole state eq. (40) 
with the second largest eigenvalue in the original seven-orbital model 
eq. (29) and the j-j coupling model eq. (30) for n=12 and £=1.0. 



dering type cannot be specified by the present calculations. 
The eigenstate with the largest eigenvalue is found to be char- 
acterized by 3g, while the eigenstate with the second largest 
eigenvalue is labelled by 2u. 

In order to examine the multipole state, in Table II, we 

(k) 

show p\ of the eigenstate with the largest eigenvalue. We 
find small difference in the values between the original seven- 
orbital and the j-j coupling models, but it can be concluded 
that the j-j coupling model works well for the description of 
the multipole state for n=l2. In the multipole state charac- 
terized by 3g, we find two significant components of rank 2 
(quadrupole) and rank 6 (tetrahexacontapole), while the rank 
4 component (hexadecapole) is negligibly small. We are in- 



clined to think that the present 3g multipole state is expressed 
by anti-bonding combination of quadrupole and tetrahexacon- 
tapole, although we cannot prove it analytically at this stage. 

In Table III, we show the components for the eigenstate 
with the second largest eigenvalue. This is the multipole char- 
acterized by 2u, which is expected to appear in filled skut- 
terudite structure due to the localized nature of electrons in 
r^ 7 . 50 ' 51) Interestingly enough, the main component is not oc- 
tupole (rank 3) as expected in the case of n=5, 50 ' 51 ' but the 
rank-7 component (octacosahectapole) becomes dominant. It 
may be concluded that pure 2u octacosahectapole occurs from 
the numerical results of both the original seven-orbital and 
the j-j coupling models. Such a high-rank multipole has been 
never observed and a way to detect it experimentally is not 
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known. However, we expect that exotic ground state includ- 
ing higher-order multipoles such as rank 6 and 7 is realized in 
Tm-based filled skutterudites. 

4. Discussion and Summary 

In this paper, we have proposed the microscopic model 
for Yb- and Tm-based compounds on the basis of the j-j 
coupling scheme. We have analyzed the impurity Anderson 
model in the j-j coupling scheme with the use of a numerical 
renormalization group technique. The results have indicated 
that the j-j coupling model works well for the microscopic 
description of the multi-/-electron state. 

Note, however, that in the present paper, we have consid- 
ered only the single conduction band in the Anderson model. 
Namely, some orbitals are assumed to be localized, but in ac- 
tual situations, / electrons in all orbitals should be, more or 
less, hybridized with conduction bands. The present Ander- 
son model is useful to select the candidates of possible mul- 
tipoles at low temperatures in an unbiased manner, although 
it is not enough to understand low-temperature properties of 
actual materials quantitatively. 

In order to consider simultaneously the formation of heavy 
electron state and the appearance of magnetism and/or super- 
conductivity, it is necessary to analyze the periodic Anderson 
model. It is not difficult to write down the Hamiltonian for the 
the multiorbital periodic Anderson model on the basis of the 
present prescription of the j-j coupling scheme for Yb and 
Tm compounds. As for magnetism, one way is to derive the 
orbital dependent RKKY interaction for the determination of 
the type of multipole ordering. Such calculations may be per- 
formed, although it is difficult to discuss the competition with 
Kondo effect. The actual analysis of the multiorbital periodic 
Anderson model in the j-j coupling scheme will be an impor- 
tant issue to be resolved in future. 

As for emergence of superconductivity, first we assume that 
heavy-electron states are formed. Then, we consider multi- 
orbital Hubbard-like model for quasi-particles with the local 
interaction in the j-j coupling scheme, which is written as 



(42) 



where is the annihilation operator for / electron at site 
i with z-component \x of j=7/2, is /-electron effective 
hopping between fi- and ^-orbitals along a direction specified 
by a connecting adjacent two sites, B is given by eq. (16), 
and I are given by eqs. (20)— (26). Here t is evaluated by 
the tight-binding approximation and it is expressed with the 
use of Slater-Koster integrals, {f fa), (//tt), (ffS), and 
(//^). 52 ' 53) The model eq. (42) will be analyzed, for instance, 
within a random phase approximation. Then, it will be neces- 
sary to proceed to the fluctuation-exchange approximation. In 
any case, we expect the emergence of exotic superconductiv- 
ity in the vicinity of ordered state. The multipole ordering can 
be also discussed in the same scheme. It is one of future di- 
rections of the research. 

We emphasize that microscopic understanding of mag- 
netism and superconductivity is one of important issues in 
the research field of strongly correlated /-electron systems. 



The microscopic research will be difficult on the basis of the 
LS coupling scheme, but it is possible with the use of stan- 
dard field-theoretical techniques if we exploit the j-j cou- 
pling model which has been shown in the present paper. We 
note that the applicability of the j-j coupling scheme is wider 
than one has naively expected from the standard textbook and 
it works even for the realistic parameter region concerning 
spin-orbit coupling and Coulomb interactions. We expect that 
the microscopic research on /-electron systems can be further 
pushed in future with the use of the j-j coupling model. 

In summary, we have proposed the prescription to con- 
struct the effective microscopic model for heavy lanthanide 
systems such as Yb and Tm compounds on the basis of the 
j-j coupling scheme. We have numerically analyzed a couple 
of Anderson models in which the local interactions at an im- 
purity site are described by using seven / orbitals and the j-j 
coupling scheme. We have found that entropy, specific heat, 
and multipole susceptibilities are well reproduced by the j- 
j coupling model. At low enough temperature, the multipole 
wave function is also well approximated by the j- j coupling 
scheme in a satisfactorily level for n=\3 and 12. 
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